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Abstract 



A master equation for the Kardar-Parisi-Zhang (KPZ) equation in 2+1 
dimensions is developed. In the fully nonlinear regime we derive the finite 
time scale of the singularity formation in terms of the characteristics of forcing. 
The exact probability density function of the one point height field is obtained 
corr esp ondingly. 
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I. INTRODUCTION 



Because of technical importance and fundamental interest, a great deal of efforts have 
been devoted to the understanding of the mechanism of thin-film growth and the kinetic 
roughening of growing surfaces in various growth techniques. Analytical and numerical 
treatments of simple growth models suggest that, quite generally, the height fluctuations 
have a self-similar character and their average correlations exhibit a dynamical scaling form. 
Numerous theoretical models have been proposed, of which the simplest nontrivial example 
is the Kardar-Parisi- Zhang equation [1] 

^-^iVhr = iyV'h + f{x,y,t), (1) 

where ^(x, t) specifies the height of the surface at point x and / is a zero-mean, sta- 
tistically homogeneous, white in time Gaussian process with covariance (/(x, t)/(x', t')) = 
2Do5(t — t')D(x — x'). Typically the spatial correlation of forcing is considered to be as a 
delta function, mimicking the short range correlation. Here we consider the spatial corre- 
lation as D(x — x') = exp(— ^^~2 ), where a is the variance of D(x — x'). When the 
variance a is much less than the system size L , i.e. cr <^ L, the model represents a short 
range character for the forcing. This same equation is believed to describe the statistics 
of directed polymers in a random medium where h has the meaning of the free energy. 
The average force on the interface is not essential and can be removed by a simple shift 
in h. Every term in the eq.(l) involves a specific physical phenomenon contributing to the 
surface evolution. The parameters ly, a and -Do ( and a) are related to surface relaxation, 
lateral growth and the noise strength, respectively. The effective coupling constant for the 
KPZ equation is given hj g = The limit 5^ oo is known as the strong coupling 

limit. Despite the intense effort in recent years, the properties of the strong coupling phase 
are rather poorly understood. Only the critical exponents of the strong-coupling regime 
{g ^ 00 OT ly ^ 0) are known in 1+1 dimensions and their values in higher dimensions as 
well as properties of the roughening transition have been known only numerically and by 
the various approximative schemes. 

Theoretical richness of the KPZ model is partly due to its close relationships with other 
areas of statistical physics. It is shown that there is a mapping between the equilibrium 
statistical mechanics of a two dimensional smectic-A liquid crystal onto the non-equilibrium 
dynamics of the (1+1)- dimensional stochastic KPZ equation 0. It has been shown |^ 
that, one can map the kinetics of the annihilation process A + B ^ with driven diffusion 
onto the (l+l)-dimensional KPZ equation. Also the KPZ equation is closely related to the 
dynamics of a Sine-Gordon chain 0, the driven-diffusion equations and directed paths 



in the random media and charge density waves ^ , dislocations in disordered sohds JTO 



formation of large-scale structure in the universe ||Tl|Jl^ , Burgers turbulence |13[ and etc. 

In this paper we attempt to clarify the short time physics of KPZ equation in the strong 
coupling limit by presenting an exact solution of the model in 2+1 dimensions. It turns 
out that eq.(l) exhibits a finite time singularity at time scale tc- Using the master equation 
method it is shown that tc can be expressed in terms of the characteristics of forcing as, 
tc = Ksc^A))^^^'^^' ^Is*^ derive the time dependence of the moments < {h — h)"^ > for 
time scales before creation of singularities. 

The main properties of the KPZ equation in the strong-coupling limit are as follow: (i) In 
the limit of ^ the unforced KPZ (Burgers) equation develops singularities for the given 
dimension. In one spatial dimension it develops sharp valleys (point shocks). In the singular 
points the height gradients are not continuous |T^. In two spatial dimensions KPZ equation 
develops three types of singularities, the first singularities are finite sharp valley lines (shock 
lines) across which the height gradients are discontinuous. The second type is the end point 
of the sharp valley lines which separates the regular points and singular region and is called a 
kurtoparabolic point. As time increases these sharp valley lines hit each other and crossing 
point of two valley lines (shock lines) produces a valley node (shock node). Generically 
kurtoparabolic points disappear at large times and only a network of sharp valley lines 
survive A complete classification of the singularities of KPZ (Burgers) equation in two 
and three dimensions, by considering the metamorphoses of singularities as time elapses, 
has been done in |[15|. (ii) For white in time and smooth in space forcing, the sharp valley 
lines (cusp lines ) are smooth curves and the singularity structures of the forced case is 



similar to unforced problem. It is shown in ||T6[ that in the stationary state the sharp valley 
lines produce a curvilinear hexagon lattice in two dimensions and, therefore, one finds a 
hexagonal tiling of singularity lines (valley lines). 

Here we shall adapt the master equation approach which enables us to investigate the 
KPZ equation in the strong coupling limit [14]. It is one of the advantages of this method 
that allows us to write all the nonlinearities due to the nonlinear term ^(V/i)^ in a closed 
form. When a is finite, the nonlinear term develops finite time singularities for the height 
derivatives in the limit {u 0). Therefore, one would distinguish between different time 
regimes before and after the sharp valley formation. Starting from a flat initial condition, 
i.e. h{x, 0) = 0,hx = u{x, 0) = and hy = v{x, 0) = 0, we aim at tracing the time evolution 
given by the KPZ equation with z/ — > 0. Only for time scales before the sharp valley 
formation {t < tc) the limit u ^ and z/ = are equivalent. For t > tc, the anomalous 
contributions of the diffusion terms to the equation of probability distribution function will 
break this equivalence [14]. Therefore, for t < tc, one can impose u = in the dynamics 



and try to find realisable statistical solutions for the problem. We show that it is possible to 
find finite time realisable solutions for the probability density function of the height field at 
least in the early steps of the evolution of the surface. We attribute the finite life time of the 
reahsable ensemble of statistical solutions to the fact that after a finite time the derivatives 
of the field h{x, t) become singular and , therefore, the diffusion term contributions are no 
more negligible. 

II. CALCULATION OF THE HEIGHT MOMENTS BEFORE THE 
SINGULARITY FORMATION 

Let us consider the Kardar-Parisi-Zhang equation in 2+1 dimensions, 

htix,y,t) - -{hl + hi) = u{h^^ + hyy) + f{x,y,t). (2) 

Defining hx{x, y, t) = u{x, y, t) and hy{x, y, t) = v{x, y, t), differentiating the above equa- 
tion with respect to x and y and neglecting the viscosity term in the limit of — > before 
the creation of singularities we get 

ut{x, y, t) = a{u{x, y, t)ux{x, y, t) + v{x, y, t)uy{x, y, t)) + /^(x, y, t) (3) 
vt(x, y, t) = a(u(x, y, t)v^(x, y, t) + v(x, y, t)vy(x, y, t)) + fy(x, y, t). (4) 

Now introducing as 

e = 

exp (^—iX{h{x, y, t) — h{t)) — iiiiu{x, y, t) — iii2v{x, y, t)^ . (5) 

enables us to express the generating function as Z{X, fii, fj,2,x,y,t) = (0). Using the KPZ 
equation and its differentiations with respect to x and y and neglecting the viscosity term, 
in the limit of — > before the creation of singularities, we get the following expression for 
the time evolution of Z{X, fj,i, //2, x, y, t) 

Zt = i-t{t)\Z - i\^{u^e) - iA^(i;'0) 

— ia^jLi{uUxQ) — iafj,i{vVxQ) — ia^2{uUyQ) — ia^2{vVyQ) 

- iX{fix, y, t)Q) - ini{fx{x, y, t)0) - iH2{fy{x, y, t)Q) (6) 

where ^{t) ^ ht ^ ^ < u'^ + >, k{x - x') = 2DoD{k - x'), k"{0, 0) = kxx{0, 0) = kyy{0, 0). 
To simplify the above expression, by taking the derivative of Z with respect to x and y, and 
using the homogeneity conditions, we obtain 



"Zx = {{-iXu - i^iUx - i^2Uy)Q) = (7) 
Zy = {{~iXv - ifiiVx - ifi2Vy)Q) = 0, (8) 

from which it can be easily derived 
d 

i- — {{-iniUx - iiJ2Uy)Q) = (u^Q) - ijJiiuUxQ) - ii22{uUyQ). (9) 
From the eqs. (|^) and (P) we have 

d 

{u^Q) - ifii{uUxQ) - ifi2{uUyQ) = -A- — (w0) = -iAZ^^^^ (10) 

O/ii 

from which we find the following expression for —iafii{uUxQ) — iaii2{uUyQ) 

— iafii{uUxQ) — iaijL2{uUyQ) = —iaXZ^-^^^ — a{uxQ)- (11) 
Similarly for —iafii{vVxQ) — iafi2{vVyQ) we find 

— iafii{vVxQ) — iafj,2{vVyQ) = —iaXZ^.^^^ — a{vyQ) (12) 



Now from eqs. ( pi] ) and (|T2|) and by using the Novikov's theorem the eq.(^) gives 

Zt = i^{t)XZ + iX-Z^^^^ + «A-Z^2M2 - i<^XZ^^^^ - iaXZ^^^^ 

- a{uxQ) - a{vyQ) - X^k{<d, 0)Z + filk"{0, 0)Z + fi2k"i0, 0)Z. (13) 

The term {uxQ) can be evaluated as follows 

1 t 
{uxQ) = — {Q)x H {{iXu + in2Vx)0) 

= -t-Z,,~^{vxe), (14) 

also for {vyQ) we have 

{vyQ) = -t-Z,,-^^{uyQ). (15) 

In the appendix it is proved that {uyQ) = {vxQ) = 0, so the equation governing the time 
evolution of Z will become 

Zt = i-f{t)XZ - iX-Zf,^^^ - iA-Z^2^2 

+ la — Z^, - la — 

- X^kiQ, 0)Z + filk"{0, 0)Z + fi2k"i0, 0)Z. (16) 



In what follows we solve this partial differential equation by using flat initial condition, 

h{x, y, 0) = u{x, y, 0) = v{x, y, 0) = 0, which means 

P{h, u, V, 0) = 6Ch)S{u)6{v), (17) 

so that the initial condition for the generating function is Z(A, /ii, /i2, 0) = 1. 
The solution of the eq. (|TB|) can be factorized as 

Z{X, /ii, ^2, t) = Fi(A, fii, t)F2{X, fi2, t) exp {-X^k{0, 0)t), (18) 

which by inserting eq.(|18D in eq.(p!6D the following equation is obtained 

tt ex 
FuF2 + F^F2t = ii{t)\F^F2 - iX-F2F,^^^^ - i\-F,F2^^^, 

+ za— F2F1 - za— F1F2 

- A'A;(0, 0)FiF2 + ^?fc"(0, 0)FiF2 + fi2k"{0, 0)FiF2. (19) 

On the other hand, 7(t) = ht = ^{u^ + v^) and ^{u'^) = ^{v'^) = -ak"{0,0)t, so 
eq.(|l9D can be separated in terms of /xi and /i2 and it can be seen easily that Fi(A, jji, t) and 
-^2(A,/U2,t) are satisfied by the following l+l-dimensional equation 

FtiX, fi, t) = -zA|f^^ + la^F^ + [fx''k"iO, 0) - ta\k"{0, 0)t]F, (20) 



which can be solved exactly with the initial condition F(A,yU,0) = 1 WM. The solution of 
eq.(|OD is 



F{fi, A, t) = {l- t&nh\J2ik^^{0, 0)aXt)) 



exp [-- ln(l - tanh.*{J2ik^^{0, 0)aXt)) 



+ -tanh~^(tanh^(Y/2zA;^^(0,0)aAt)) - X^k{0,0)t 



- — In 



1 , 1 -tanh(y'2zA;,,(0,0)aAt), 



16 1 + tanh(./2iA;^^(0,0)aAt) 



-^z/x^/ ^'^""^°'°^ tanh(722fc..(0,0)aAt)]. (21) 
2 V aX ^ 

So, the solution of eq.(|l9D is 

Z(A, /xi, /i2, t) = F(A, /xi, t)F(A, /i2, t) exp {-X^k{0, 0)t) 



(1 - tanh2(y'2iA;^^(0, 0)aXt)y exp[-- ln(l - tanh^(y^2iA;^^(0, 0)aXt)) 



+ -tanh ^(tanh^ 



2i/cxx(0,0)aAt)) - A^/c(0,0)t 



1 - tanh(./2iA;^^(0,0)aAt), 



tanh{J2ik:,^{0,Q)aXt) 



1 2^ 2 /2zM0,0)^ 



tanh(A/2iA;^^(0,0)aAt)]. 



(22) 



Using the generating function, the probabihty distribution function (PDF) of height 
fluctuation can be derived by inverse Fourier transformation as 



P{h, u, V, t) = 
dX dfii dii2 



I 



exp (iXh + ijjLiu + ifi2V^Z(X, /xi, //2, t)- 



(23) 



27r 27r 27r 

By expanding the solution of the generating function in powers of A, all the {{h — h)"") mo- 
ments can be derived. For instance, the first five moments before the sharp valley formation 
are 

2/3 r 

t 



ak" {0,0) ) 
48 / P(0,0) 



45 \ak" {0,0)^ 



ak" (0,0) J 

( k\0,0) ' 
\ak"{0,0)^ 



35 H* 



5/3 



1216 t 
945 ^t*' 



+ 



64 t 



(24) 



where — (^sfr^^o^)^^^- An important information content of the derived exact form is 
that they determine the time scale of the singularity formations. One should first check the 
realisability condition, i.e. P{h — h,t) > 0. In fact the moment relations above indicate that 
different even order moments will get negative in some distinct characteristic time scales. A 
closer inspection into the even moment relations reveals that the higher the moments are, the 
smaller their characteristic time scales become. So the rate of decreasing tends to tc — |t* for 
very large even moments asymptotically, where U = (^^^^tS^^)^''^ = ( g^^^^ )^/^(T^. Therefore, 
it is concluded that at this time the right tail of the probabihty distribution function (PDF) 
starts to become negative, which is the reminiscent of singularity creation. Physically, 
the connection between the negativity of the right tail and the singularity formation time 
scale is related to a multi-valued solution of KPZ (Burgers) equation for t > tc- It should 



be notified that eq.(16) preserves the normahsabihty, i.e. Z(0,0,0,t) = 1, equivalent to 
J-^ P{h — h,u,v;t)d{h — h)dudv = 1 for every time t. So the PDF of h — h and its 
derivatives are always normalizable to unity. In the u = limit, for time scales greater 
than tc , the height field will become multi- valued in the valleys (the shock region), which 
is related to the P{h — h) left tail . The multiplicity of the height field will increase the 
probability measure in the PDF. Therefore, to compensate the exceeded measure related 
to the multi-valued solutions the right tail of the PDF tails should become negative. The 
singularities in the limit z/ — can be constructed from multi-valued solutions of the KPZ 



equation with = by Maxwell cutting rule |T^, which makes the discontinuity in the 
derivative of the height filed. When t > tc, the contribution of the relaxation term in the 
limit of vanishing diffusion coefficient should also be considered in order to find a realisable 
probability density function of height field. In other words, the disregardance of the PDF 
equation's diffusion term is only valid up to the time scales in which the singularities have 
not been developed yet. 

Taking into account a > and k"{0,0) < 0, the odd order moments will be positive 
in time scales t < and the probability density function P{h — h,t) is positively skewed. 
Therefore, the probability distribution functions of height field will have a non zero skewness 
through its time evolution at least before the singularity formations time scale. The inverse 
Fourier transform of the generating function has been performed numerically to obtain the 
PDF P{h — h) form. To demonstrate the time scale of the singularity formation the PDF 
time evolution has been sketched numerically, Fig(l). As the system evolves in time, the 
formation of the first singularities leads to the right tail negativity in the PDF. 

To summarise, we analyse the fully nonlinear KPZ equation in 2-1-1 dimensions forced 
with a Gaussian noise which is white in time and short range correlated in space. In the 
non- stationary regime when the sharp valley structures are not yet developed we find an 
exact form for the generating function of the joint fiuctuations of height and height gradient. 
We determine the time scale of the sharp valley formation and the exact functional form of 
the time dependence in the height difference moments at any given order. In this paper, 
we have used the initial condition h{x, 0) = and u{x, 0) = 0. This method enables us to 
determine the finite time scale of singularity formation in terms of statistical properties of 
the initial conditions. This calculation will be dealt with elsewhere. At this level we would 
remark on the stationary state of the KPZ equation in the strong coupling limit. Using 
the numerical results recently, it is shown that convergence to the statistical steady state is 
reached after a few turnover times [16]. Therefore, in the stationary state the singularities 
will be fully developed and, therefore, at large t or t — > oo we should take into account the 
relaxation contribution in the limit z/ ^ 0. Contribution of the relaxation term in the rhs 
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FIG. 1. In the LHS graph the time evolution of PDF of h — h before singularity formation at 

and is numerically obtained. RHS graph shows the right tails of the PDF of h — h for ^t* 
and corresponding to before and after singularities formation which are numerically calculated. 



of the eq.(6) is, 



lim{-zzy/ii( V^uO) - zzy/i2( V\e)}. (25) 



In we have calculated the finite contribution of this term in 1+1 dimensions and 
showed that it plays an important role in statistical theory of KPZ equation in the station- 
ary state. It is also shown that all of the amplitudes of the moments < {h — /i)" > can be 
expressed in term of the relaxation term in the limit u —>■ 0. In 2+1 dimension calculation 
of finite contribution of the relaxation term is more complex ( because of complex structure 
of the singularities) and it is left for a future work. We believe that the analysis followed in 
that paper is also suitable for the zero temperature limit in the problem of directed polymer 
in the random potential with short range correlations [10]. In the same direction the present 
method is applicable to the strong coupling regime of KPZ equation in higher dimensions 
{d > 2) which is definitely an important step. 
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III. APPENDIX 

In this appendix we present another way of deriving the height moments by introducing 
the second spatial derivatives of height field. Defining 6(A, fii, /i2, Vi^V2, Vsi Vj ^) &s 

e(A, /ii, /i2, r^i, ??2, X, y, t) = exp {-iX{h{x, y, t) - h{t)) - ifiiu{x, y, t) - ifi2v{x, y, t) 

-ir]iw{x, y, t) - ir]2s{x, y, t) - i'n3q{x, y, t)}, (26) 

the generating function is defined as 

Z{X,fii,^i2,rii,r]2,r]3,x,y,t) = (6), (27) 

where u{x,y,t) = h^{x,y,t),v{x,y,t) = hy{x,y,t),w{x,y,t) = hccx{x,y,t), s{x,y,t) = 
hr,y{x,y,t) and q{x,y,t) = hyy{x,y,t). 

The evolution of h{x, y, t), u{x, y, t), v{x, y, t),w{x, y, t), s{x, y, t), q{x, y, t) is given by the 
following equations 

h, = ^{u' + v') + f{x,y,t) (28) 
ut = a{uw + vs) + fa:{x,y,t) (29) 



Vt = a{us + vq) + fy{x,y,t) (30) 

Wt = a{w'^ + + uwa, + vsa,) + fxx{x, y, t) (31) 

St = a{ws + qs + vSyUWy) + fxy{x, y, t) (32) 

qt = q;(s^ + + uSy + vqy) + fyy. (33) 

It follows from the above equations that the generating function Z is the solution of the 
following equation 

CK Oi 

Zt = i^{t)\Z - i\-Zn^n^ - iA-Z^2^2 + + ^«^m22/ 

- iaZjj-^ — iaZri^ + iarjiZ,i-^,i-^ + iarjiZ^^r,2 

+ iar]2Zr,^r,2 + i(^V2Zn2m + ^«^3^r,2»72 + ^"^3^r,3ry3 

- ^A(/(.'z:, 2/, t)0) - ifiiifxix, y, t)Q) - iiJ,2{fy{x, y, t)Q) 

- iVi{fxx{x,y,t)Q) -ir]2{fxy{x,y,t)Q) - ir]3{fyy{x,y,t)Q), (34) 

in which j{t) is defined as 7(t) = /it and the following identities have been used 

r]i{wxQ) + m{sxQ) + VsiQx^) = iZx - i\Z^^ - i/iiZ^^ - i//2^„2 (35) 
r]i{wyQ) + r]2{syQ) + Vsily^) = ^Zy - iXZ^,^ - ifiiZ^^^ - ifi2Zrj3. (36) 

Now, using Novikov's theorem we find 

{fix, y, t)e) = -iXk{0, 0)Z - ir]ikxx{0, 0)Z - imk^^iO, 0)Z (37) 

{fx{x,y,t)e)^ifXik,,{0,0)Z (38) 

{fy{x,y,t)e)^ifX2kxx{0,0)Z (39) 

{fxx{x, y, t)e) = -iXka:x{0: 0)Z - ir]ikxxxx{^, 0)Z - ir]3k^xxx{0, 0)Z (40) 

{fxy{x, y, t)Q) = -ir]2kxxxx{0, 0)Z (41) 

{fyy{x, y, t)Q) = -iXk^xiO, 0)Z - ir]ikxxxx{0, 0)Z - irjskxxxxiO, 0)Z, (42) 

where k{x -x',y- y') = 2D^D{x -x',y- y'), k{0, 0) = kxx{0, 0) = kyy{0, 0) ^ 
and kx{0, 0) = ky{0, 0) = 0. So we have 

= i^{t)XZ - iA^Z^i^i - iA^Z^2/U2 + ioiZi^^x + iaZ^^y 

- iaZjj^ — iaZjj^ + iarjiZ^^^^ + iai^iZ^^^ 

+ ia'q2Zrfir)2 + io(r]2Zr,^r,3 + ioLr}2,Zr,2r,2 + ioir]sZ,j^,j^ 

- X'^kiO, 0)Z + ii^l + 14- 2Xr]i - 2Xr]s)kxx{0, 0)Z 

- ivl + Vl + Vl + '2r]ir]3)kxxxx{0, 0)Z. (43) 



Assuming statistical homogeneity {Z^ = 0,Zy = 0) and defining P{h,u,v,w,s,q,t) as 
the joint probabihty density function of h, u, v, w, s and q, one can construct the PDF as the 
Fourier transform of the generating function Z with respect to A, /ii, 112, Vij V2, 

~ f dX dni dfi2 drji drj2 drj3 

P{h,u,v,w,s,q,t) ^ / — ^ ^ ^ ^ ^ { 



27r 27r 27r 27r 27r 27r 

exp (j,Xh + i/J-iu + iiJ,2V + ir]iw + ^772^ + i?73g)z(A, //i, /X2, r/i, 772, ?]3, 0}- (44) 

From the eqs. (43) and (44) the equation governing the evolution of P(h, u, v, w, s, g, t) 
can be derived, which is 

= 7(i)P;, + 17(«^ + 'y^)-Fk - 4awP - AaqP 

— avj^Pw — a(fPq — as^Pw — as^Pq — awsPg — aqsPg 

+ A;(0, 0)P;,^ + 2A;,,(0, 0)P^^ + 2A;,,(0, 0)P^^ - A;,,(0, 0)P„„ - A;,,(0, 0)P,, 

From the eq.(45), it is easy to see that the moments {h^°u^^v'^'^w'^^s'^^q^^) satisfy the 
following equation 



dt 



ario 



+715(715 - l)h,,,{0,0){h'''u'''vV's''^q'''-'') + 2nsn5h^^^{0,^^ 

+ri3(rz3 - l)A;^^^^(0,0)(^"°'u"i'y"2w"3-^s"*g"«) +?^4(?^4 - 1)/Cxxxx(0, 0)(^"°M"i'f;"^'u;"3s"*-2g"'^) 
+an5(/i""«"^'y"^w"3s"^+^g"^-i) + tiqK - 1)A;(0, 0)(/i"°-2«"^'y"2w"3s''^g"«). (46) 

By substituting different sort of values for no,ni,n2,n3,n4^, 715 we can find some coupled 
differential equations for different moments. For example, if we take no — rii — ■ ■ ■ — — 
we have 

<w> + <q>=0 (47) 

or 

< V • u >= (48) 

where this is the same as the statistical hemogeneity condition. For uq = l,ni = n2 = = 
7i4 = 72.5 = we find 

<hw > + <hq>^ - <u^ > - <v'^ > (49) 
and for no — 0,ni — l,n2 — n^ — n^ — n^ — we find 

< uw > + < uq >^ 0, (50) 
and if we assume the statistical hemogeneity we have 

< UW > = < UUx >= 2 ^ ^'^ >x= 0. (51) 

So < UW >— and then < uq >— 0, also for 710 = ?^! = 0, n2 — l,ns — n^ — n^ — we find 
< vq >— 0, < vw >— 0. 

For no — 0,ni — 2,n2 — n^ — n^ — n^ — we find 

^ < i(2 -a < u'^w > -CK < u^q > -2kxx{0, 0) (52) 



and also we have < u^w >— ^ < >x= by statistical hemogeneity and also 

< u^q >—< u^v >y —2 < uvs >, < u^v >y= 0, so < u^q >— —2 < uvs > so that we have 

^ <u^ >^2a< uvs > -2k^^{0, 0). (53) 

The corresponding differential equation for < uvs > is 

d , . 

— < uvs >= 0. (54) 
at 

If we assume that a,t t — the surface is flat so all the moments a,t t — are zero, so that 

< uvs >— and then we find 

Ku"^ >^ -2k^^{0,0)t. (55) 

Similar calculations give 

< v'' >= -2k,,{0, 0)t. (56) 

By focusing on the differential equation of the PDF it is deduced that this equation is 
invariant under u —u and v — > —v which is a consequence of the inversion and reflection 
symmetry of the KPZ equation. So it will result in the vanishing of the u and v odd moments 

Using eq.(46) the even moments of u and v can be calculated. For example, for {u^) we have 

^U) = -a{u'w) - a{u\) - 12/c,,(0)(«2). (58) 
at 

On the other hand (u^w) can be written as 

= (59) 
which is zero by homogeneity, also for {u'^q) we have 

{u\) = {u^v)y - A{u^vs) = -A{u^vs) (60) 
in which {u%)y is zero by homogeneity. The differential equation for {u^vs) is 

j^{u'vs)^Q, (61) 
so {u^vs) = 0. Therefore, it is obtained 

(«^) = mL(o, o)^^ (62) 



and we have the same for {v'^). Also it can be found that {u^v"^) = 4A;^^(0, 0)t^. By continuing 
the above method all the moments {u^'v™-) can be deduced. 

On the other hand, the above calculations show that all the mixed moments {u^v"'s) are 
zero. In the following some of the moments of u and v and their combinations have been 
given 

{u') ^ {v') ^ -120kl{0,0)t' (63) 
(«V) = («V) = -24/.L(0,0)i^ (64) 

{u^) = {v^) = lQ80kU0,0)t'' (65) 
(«V) = (mV) = 240A:^,(0, 0)t^ (66) 
{u%^) = 144A;^^(0, 0)^^ etc. (67) 

We aim to calculate all the momentes oih — h — h. Putting hq — 2, rii — n2 — • • • — 
in the eq.(46), we find 

A(/i2^ ^ -a{hu^) - a{hv^) - aiJi'w) - aiJi'q) + 2k{0,0). (68) 

We know that {ft^w) — {h^u)x — 2{hu^) and because of the homogeneity {h^u)x — 0, so we 
have 

iji^w) = -2{hu^), (69) 

and also 

{h\) = -2{hv'), (70) 

then 

^(h^) = aUhu^) + (hv^)) + 2k(0, 0). (71) 
at 

As it can be seen from eq.(71) we need the moments (hu^) and (hv'^) to find {Jt^) , which 
the related differential equations are 

^{hu^) = -j{t){u^) - ^{u") - ^{u%^) - a{hu^w) - a{hu'q) (72) 

f/hv') = -^{t){v') - ^{v') - ^{u\') - aChv'w) - a{hv\). (73) 

By using of the statistical homogeneity the last two terms of the above equations can be 
converted as 



(hu^q) = —{v^v^) — 2{huvs) 

(hv^w) = -{u^v^) - 2{huvs). (74) 



The above relations result in 



.2„.2 



+ a{u^v^) + Aa{huvs). (75) 

Obtaining the differential equation for (huvs) we get 

^(huvs) = 0, (76) 

which results in (huvs) — 0. Also it is easy to see that all the moments {h^u^v^s) are 
zero. This fantastic result helps us to find all the (k^) moments. Now by substituting 
7(i) ^ht^ f ((li^ + v^)) = -'^O(k^x{0, 0)t , {u'^) , {v*) and {uH'^) in eq.(75) we obtain 



which finally gives 



{hu') + Chv') ^ --akl,{0,0)t', (77) 



{h') = —a^klMt' + 2kMt. (78) 



Now we begin to calculate the moment (h^)- Inserting riQ — 3,ni — n2 — • ■ • — — 
in eq.(46) we get 

-{h^) = -3j{t){h^) - ^a{{h^u^) + {h\^)) - a{h?w) - a{h?q), (79) 

and again by using statisiacal homogeneity it can be shown that 

(h^'w) = -Z(h^v?) (80) 

(h^q) = -?>(h^v''). (81) 

So we have 

= + ^«((/^v) + (/.V)). (82) 

To calculate Qv") we need (^) , {J^v} + Jv^v^)). (fi^) has been calculated above, so we will 
obtain iji^v? + fv^v^) using the corresponding differential equation as follows, 



^{{h\^) + (h\^) = -2j{t){{hu' + hv^)) - a{{hu^ + hv^)) - 2a{hu\^) - a{{hWw + h^v^w)) 
- a{{h^u^q + h%\)) + 2k{0, 0){{u^ + v^)) - 4A;^^(0, 0){h^). (83) 

As before, we easily get 











(hWq) = 










-2{hu\'') 


- 2{h'^uvs) 






-2{huW) 


-2{h^uvs). 


(84) 



As discussed before, it is easy to show that the moment {hi^uvs) is zero. To prove this 
we write the corresponding differential equation 

^iji^uvs) = -a{{hu^vs) + {huv^s)), (85) 
at 

and again by trying to write the differential eqautions for (hu^vs) and (huv^s) we obtain 

^iWvs) = 
at 

^{huvh) = (86) 
at 

which results in (hu^vs) = (huv^s) = 0, therefore (h'^uvs) = 0. Now (hu^) , {hv"^) and 
{hu^v"^) should be found. The relating differential equation for {hu'^ + hv"^) is 

f^iihu') + {hv') = -j{t){{u' + v'))- ^{{u' + v')) - |((^V + uV)) 

- a{{hu^w + hv^w)) - a{{hu^q + hv^q)) - 12A;^^(0, 0){{hu^ + hv^)). (87) 

As before the following identities are held 

{hu^w) = --(u^) 
5 

(hu^q) = -(^^V) 

(hv'q) = -^{v') 

(hv^w) = -(mV), (88) 



so we have 



- 12A;^^(0, Q){{hu^ + hv^)). (89) 



Substituting the expressions for 7(t) ,{u^) , (v^) , (u^) , (y^) , {u^v'^), {v?v^) and {{hu'^ + hv'^)) 
we find 

{hu'^ + hv^) = S2ak^^(0, 0)H\ (90) 
The corresponding difi^erential equation for {hu'^v^) is 

|((^uV) = -7W((«V)) - |((^.V + «V)) 

- a{{hu\^w + hu^^q)) - 2k^^{0, 0){{hu^ + hv^)). (91) 

As before the following identities are held 

(hu'v'q) = -^(kV), (92) 
then by inserting the known moments we find 

{hu'v') ^ ^aklMt\ (93) 

The above results get 

+ {h'^v^) = -^a^kl^iO, 0)t' - 8ki0, 0)A;,,(0, 0)^^ (94) 



which finally results in 



48 

{h') ^ --aX^Mt^ (95) 



By continuing the above proccedure all the (/i") moments can be derived. Some of these 
moments are listed bellow 

44 

{h') = -^c^XM 0)^' - Sa'k{0, 0)kl{0, 0)t' + 12k\0, 0)t' (96) 

ih') = -^«'^L(0, 0)t'" - fa'k{0, 0)klM 0)f, (97) 

which are the same as the results that we derived by expanding the generating function. 
Results of this appendix show that any moment containing the first power of s vanishes. 
Indeed, one can prove the following identity 

^^^-i(xh+nM2v)s^ = 0. (98) 



By expanding the exponential in the above expression one finds 



Jn+m+p) \n . m . P 

^ ^ nj^,p ri\m\p\ ^ ^ ^ ^ 



Now setting — — and 77.4 = 1 in eq.(46) we get the following equation for 

-p{p-l)K,{0,0){h''u"'vP-h). (100) 

Starting from eqs.(54),(61) and (76) all the (sh'^u'^v^) moments can be evaluated and 
it can be shown that, with flat initial condition, they are equal to zero. So we conclude 
{sQ) = {uyQ) = {v,Q) = 0. 
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